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' Abstract. Increasing observational evidence supports a picture for a 

^ . close link among AGN phenomenum, star-formation processes and galaxy 

formation. Since the physical phenomena related to the onset of both 
AGN and star formation are very likely characterized by a strong far 
I infrared (FIR) emission it becomes mandatory to investigate the FIR 

energy domain. However, because of their faintness in the FIR energy 
range, very little is known about the FIR properties of type 1 AGN, and 
type 2 AGN, which are supposed to be missed in classical searches, still 
^T) I lack definition and samples. 

Our aim is to investigate the FIR properties of optically selected type 1 
AGN, compare them to Ultraluminous Infrared Galaxies (ULIRGs) and 
derive some general characteristics of the population. Expectations for 
[ future surveys are also presented. 
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For historical reasons and selection criteria, a classical type 1 AGN is that show- 
ing a spectral energy distribution (SED) with most of its energy in the UV and 
only 10 - 20 % in the infrared. From 10 fim to 1000 A the agreement among 
■ different composite spectra, built on observations from different databases, is 

\ excellent, while in the FIR the various studies disagree. The FIR SED differs 

not only in shape but also the FIR luminosity differs by up two order of mag- 
nitudes. This difference is due to inhomogeneities of AGN samples, covering 
different luminosity ranges and objects selected according to different criteria. 
However a common result of these studies is that both AGN and starburst (SB) 
dominated objects have a common FIR-emitting star-formation component from 
which is not easy to disentangle pure AGN emission. Quite a number of obser- 
vational evidence supports the above claim: (a) their 60 fim flux, radio and 
CO emissions are well correlated and AGN and star-forming galaxies show the 
same properties; (b) their SEDs at A > 60 /im have the same shape as those of 
star-forming regions; (c) the FIR Luminosity Function of AGN coincides with 
that of SBs in shape and differ only in amplitude: ^agn/^^b ~ 1/^5; (d) the 
evolution rate of the two classes is comparable and shows the same behaviour 
L{z) cx (1 -|- z)^"^ (e) the FIR component represents a distinct component cor- 
relating only weakly with the true AGN emission and very likely related to a 
concurrent SB (Lawrence, 2001; Rowan-Robinson, 2001). 
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All this supports the idea that all AGN are accompanied by a SB but one SB 
out of 25 is accompanied by an AGN. There are indeed evidences for a coexis- 
tence of central AGN and circumnuclear star formation in a significant fraction 
of ULIRGs (Genzel et al., 1998) but AGN are difficult to be observationally 
proved mainly because of the high extinction of the nuclear source. Other dis- 
tinct direct and indirect observational probes indicate that an unknown fraction 
of AGN is absorbed, the so-called type 2, and the ratio between type 2 and type 
1 may be as high as 10 of which 20 percent in FIR selected objects. For in- 
stance, the comoving mass density of Massive Dark Objects, which are probably 
supermassive Black Holes in galaxy centers, are a factor of 2 or larger than that 
explained by accretion onto optically-luminous AGN. This implies that much 
accretion is obscured by dust (Haehnelt, Natarajan and Rees, 1998) a popula- 
tion which is missing in optically selected samples. Other evidences come from 
X-ray studies and are extensively discussed in this volume. 
Dust enshrouded type 1 AGN exist at high and low redshifts and constitue a 
subset of SCUBA galaxies, which in turn share common properties with local 
ULIRGs. Nearly all nearby ULIRGs appear to be late-stage mergers, then if one 
out 10 SB has quasar-like activity is then possible to envisage the evolutionary 
path starting from a merging system going through a starburst and ending to a 
optical-phase QSO (Sanders et al., 1989). 



2. Present work 

The aim of the present work is to use FIR observations, gathered by the Infrared 
Space Observatory (ISO), of optical AGN to infer any evolutionary link from 
'pure' AGN objects to eventual transition objects linking AGN to ULIRGs. 
This aim was pursued by comparing FIR and optical properties in a complete 
subsample of optically selected bright (15 < B < 17) quasars. The sample 
completeness ensures a better control of observational biases and provides reli- 
able information on the properties of the underlying population. The objects 
studied in this work can be then fully considered as representative of the entire 
population. The optical observations consist of multicolour (U,B,V,R,I) broad- 
band photometry used primarily to selected quasars candidates, which were 
then later spectroscopically confirmed. ISO FIR photometry has been taken at 
7,12,25,60,100 and 160 fim. All the details of these observations and the sample 
properties can be found in Andreani et al. (2001). 

3. A colour-colour diagramme for QSOs 

Figure ||(A) reports a(60/xm,MIR) versus a(FIR, 60/xm). By MIR we indi- 
cate the mid-IR bands 7, 12 and 25 nm and by FIR the far-IR bands at 100 
and 160 ;um. a are the spectral indices, i.e. for a /jy oc z^" SED: 0(1^1,^2) = 
log(/(i/i)//(z^2)) — log(i^i/i^2) and are therefore simply related to colours. Our 
results (see Andreani et., 2001) indicate that a(FIR, 60/irn-) is a fair estimate of 

where T„d is the temperature of a warm dust component. 
The diagram shows a trend: the relative emission at MIR (shorter) wavelengths 
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with respect to the 60//m emission decreases as T^d increases. A possible inter- 
pretation is the increasing dominance of SB component with respect to AGN in 
objects with lower values of T„d and lower values of the ratio F(60/im) /F(MIR). 
The colour-colour diagram was indeed widely used in the past as a tool to detect 
and discriminate between different types of activity in the nuclear and circum- 
nuclear regions of galaxies (see e.g., Canalizo & Stockton (2001) and references 
therein). Different kind of objects (QSOs/Seyferts, starbursts and ULIRGS) oc- 
cupy distinct areas of such a diagram: where FIR fluxes are dominated by dust 
reradiation objects locate in the lower right corner (in our convention at lower 
T„d and lower Q(60^m,MIR)) while those with strong non-thermal emission, 
such as optically selected QSOs, are in the upper left corner. This means that 
our sample spans a wide range of properties from those 'AGN pure' to those 
where a concurrent SB dominates. 



4. The type 1 SED 

By making the rather coarse assumption that optical quasars can thus be mod- 
eled as a homogeneous population, observations at fixed frequencies of targets at 
different redshifts can be then used to build the entire SED for this population. 
This analysis exploits the whole spectral information from the optical to the 
FIR and helps outlining possible general features of the emitting mechanisms 
and the physics of the QSO environment. 

The spectrum is built by dividing the available wavelength range in bins. In 
each bin an average of the fluxes at that wavelength is computed. The errorbars 
shown are those related to the averages but take into account also photometric 
uncertainties. Figure |l|(B) reports the average spectrum in the QSO restframe. 
A well defined IR component peaking at 10-30 /im and dropping steeply above 
100 fim is evident. Although from the characteristic shape it is quite straight- 
forward to infer a thermal origin for this component it can be due either to a SB 
emission in the host galaxy and/or to a dusty torus around the central source. 
When comparing these data with models (Granato & Danese, 1994) it is possible 
to bound the available parameters on the dust distribution. It turns out that 
the dust distribution around the central source should be compact because of 
the sharp cut-off at long wavelengths of the FIR bump. A larger bump would 
imply larger temperature distribution and very likely larger spatial dimension 
of the dust distribution. Furthermore, the data are better in agreement with 
models predicting the strongest dust self-absorption around 10-30 ^m. 



5. Expected counts of Type 1 AGN 

Based on these results we show here the predictions of simple models on the 
expected numbers of type 1 QSOs in deep future surveys such those foreseen 
with SIRT-F and HERSCHEL Satellites. 

As well known number counts of extragalactic objects represent one of the fun- 
damental tool to investigate galaxy evolution, since they are related to the Lu- 
minosity Function. The differential counts, dN{S)/dS, can be expressed as a 
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Figure 1. Left (A): The ratio between 60//m flux and MIR fluxes is 

shown as a function of the ratio between FIR and 60 fim fluxes. Aster- 
isks correspond to the 25 /im data, filled squares to 12 /im, open circles 
to 7 /xm ones. Open squares correspond to QSOs at z < 2 reported in 
Andreani, Franceschini & Granato (1999). Distinct areas can be iden- 
tified from AGN dominated (in the upper left corner) to SB dominated 
objects (in the lower right corner), according to the expected dominant 
emitting mechanism in the object SED from non-thermal to thermal. 
Right (B): The average spectrum in the QSO restframe. The spectrum 
is normalized at Abiue = 0.44 /xm. Filled circles correspond to the op- 
tical and near-IR data, open triangles to the ISO measurements. The 
lines going through the data show the predictions of models by Granato 
and Danese (1994). The solid line corresponds to a face-on compact 
configuration of the torus with absorption at 30/Ltm of r=30; the dot- 
ted and dashed lines to a 45° inclined torus with larger dimension and 
T=30 and r=60, respectively. 
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simple integral of the epoch-dependent luminosity function, ^{L,z): 

dN /-^^ dV dL{S,z) ..^ , . , , . 

d5=i sources/sr (1) 

where Zu and zi being the effective upper and lower limits of the redshift dis- 
tribution and ^ the differential comoving volume. The integrated counts are 
easily found by integrating eq.(l) at different flux limits, S. 
As ^{L,z) we take the B-band luminosity function that computed by Boyle et 
al. (2001) and as its redshift evolution that described in Bianchi, Cristiani &: 
Tae-Sum (2001). The latter has a down-turn at z = 3 and decays exponentially 
until z = 10. The cosmological parameters used to compute the differential 
volume are Ho=50 km/s/Mpc, 17 = 1, A = 0. We compute k-corrections, i.e. 
dL{S, z)/dS , according to the composite spectrum of Figure 0(B) and plot in 
Figure ^ the expected counts in different observing bands: 5, 8, 24 and 60 /im. 
i.e. those corresponding to the filter sets of IRAC and MIPS on board of SIRT-F 
and PACS on board of Herschel. It is easily shown that for a deep survey in 
a sky area of O.Sdeg^ at a depth of 0.4 mJy at 24;um and of 34/iJy at 8/im we 
expect to see with SIRT-F a few tens and a few hundreds objects, respectively. 
Note that the numbers presented here refer to an ideal survey not affected by 
observational uncertainties (confusion noise, photometric errors, ...). More re- 
alistic predictions taking into account also the selection of candidate QSOs will 
be the subject of a forthcoming paper. 



Table 1. Expected counts of type-1 AGN in future IR surveys 



A (/um) 


Smin(^Jy) 


sky area {deg^) 


^ of objects 


Instrument 


4.5 


5 


0.3 


1000 


IRAC/SIRT-F 


8.0 


34 


0.3 


350 


IRAC/SIRT-F 


24 


370 


0.3 


50 


MIPS/SIRT-F 


75 


5000 


1.0 


22 


PACS/HSO 
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Figure 2. Expected number counts of type 1 AGN computed on the 
basis of the composite SED of the type 1 AGN of the sample presented 
in this work. The different wavelengths correspond to the bands of 
IRAC and MIPS instruments on board of SIRT-F and PACS on board 
of Herschel. 12/L(m counts are only shown for comparison. 



